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1.0 INTRODUCTION AND SUMMARY 


The primary purpose of the World Surface Map is to provide the LAMMR 
subsystem with world surface type classifications that are used to set up 
LAMMR Level II process control. This datr base will be occessed solely 
by the LAMMR subsystem. The SCATT and ALT subsystems will access the 
data base Indirectly through the T^^ (Brightness Temperature) Data Bank, 
where the surface types have been updated from a priori to current clas- 
sification, and where the surface types have already been organized on an 
orbital subtrack basis. 

The single most important factor in the design of the World Surface 
Maps is the ease of access to the information while the complexity of 
generating these maps is of lesser importance because their generation is 
a one-time, off-line process. The World Surface Map provides storage of 
information with a resolution of 7 km necessary to set flags concerning 
the earth's features with a different set of maps for each month of the 
year. 

• Land; 

• First year ice; 

• Multi-year ice; 

• RFI at 6.6 GHz, 5.1 GHz and 4.3 GHz, and 

• Ocean. 

Properties of the World Surface data base of importance are; 

• a "telescoping" scheme to minimize total data storage, 

• geographic contiguity, 

• access efficiency, and 

• uniform distribution of cell centers. 

Several earth-oriented data base structures have been conceptual- 
ized, in numerous studies (References 1, 2, and 3), including the spher- 
ical icosahedron and the spherical cube. These various methods are con- 
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sidered for adaptation to the LAHMR World Surface Map data base struc- 
ture. Indexing of these data bases, transformation between geographic 
coordinates and computer storage Indexing, data accessing, and retrieval 
have been examined. 

Results of this study show that a latitude-longitude Indexing system 
Is most suitable for the World Surface Map. This Indexing system elimin- 
ates both the necessity to project a geometrical structure onto a circum- 
scribing sphere, and the need to explicitly compute Indices for every 
cell. Thus, retrieval Is simple and rapid. The LAMMR data Is In Instru- 
ment coordinates, with latitude-longitude addressing tags accompanying 
the data. 

The main features of the World Surface Maps are summarized as follows: 

* The earth's surface Is divided Into a number of 
geographical ly-contiguous regions, each associated with a 
data base record. 

* The addressing scheme Is binary sequential. This permits 
"telescoping" whereby regions of diminishing size are 
addressed by adding bits to Indices of Increasing length. 

* The cell addresses are binary coded latitude-longitude. 

Conversion from one coordinate system to another Is not 
necessary. 

To achieve approximately 7 km gridding Intervals a total storage of 
about 16 million cells Is required for each of the latitude-longitude 
based World Surface Maps. The earth's surface Is partitioned Into 2N^ 
squares with N equidistant partitions on the latitude axis and 2N 
equidistant partitions on the longitude axis. Each square Is considered 
a logical record consisting of M x M cells. Each cell has a gridding 
Interval of 7 km along the latitude axis. The gridding Interval along 


1-2 


Also Technoukkal SvsTEm, Inc 


the Equitor Is also 7 km, but for any other parallel of latitude, the 
gridding Interval Is less than 7 km. Results of tH' study reported In 
Section 2 show that N ■ 256 and M ■ 12 Is a good choice and N can be 
further broken down Into telescoping levels to reduce the amount of stor- 
age required. 

Seventeen bits are required for Indexing one of the 2 x 256^ ■ 
131072 regions. Another 8 bits are used for addressing cells within the 
region. Coordinate transformations are not necessary because both the 
geographic coordinates and the data base Indices are the same. However, 
there Is a requirement to encode latitude-longitude coordinates Into 
binary addresses. 

By Incorporating telescoping Into the data base retrieving process, 
the amount of storage required Is reduced. A rough estimate of 25% 
reduction In storage Is possible; the data base Is thereby reduced to 12 
million units per map. Four bits per flag are needed for storage of the 
Information. The total storage for 12 such maps Is 576 million bits. 

It should be noted here that a geographically-referenced World 
Brightness Map may be required for data comparisons between different 
spacecrafts or Instruments as well as data gathering for research pur- 
poses. The World Brightness Map requires cells of equal area and minimal 
distortion. This requirement Is most likely to be fulfilled by spherical 
Icosahedron structure, the details winch are discussed In Section 3. 

1.1 Future Study and Recommendations 

The World Surface Map needed to provide the surface type classifica- 
tion flags has been organized In area blocks within area blocks In a tel- 
escopic fashion. The purpose of this Is to provide an efficient means of 
access to this data bank and to reduce the total data volume required. 
Total data volume requirements and efficient access schemes are Important 
considerations In the design of the LAMMR subsystem algorithms. Addl- 
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tional study Is required to determine whether the selected method or some 
other such method Is most applicable. At the same time, with the advanc- 
ing computational and storage capabilities of ground systems, such as 
will be lirolemented for the PPF (Primary Processing Facility), data vol- 
ume might play a loss significant role than anticipated. As a result, 
the telescopic approach to the World Surface Map may turn out to be not 
worth the overhead associated with It. This question can only be firmly 
resolved after the hardware system of the PPF has been defined. 
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2.0 DATA BASE BASED ON LATITUDE»LONG!TUDE COORDINATES 


This section presents results of the letitude-longitude coordinates 
Earth data base study. The gridding system of this data base Is based on 
the latitude-longitude geographical coordinate reference system on the 
spherical Earth. In this coordinate, the latitude and longitude are 
treated as the X and Y axis respectively by mapping the Earth on a 
plane. The Earth as a plane reference surface with latitude-longitude 
gridding system Is shown In Figure 2-1. The gridding Intervals along the 
X-axis, the meridian, are equal. At first sight this may seem to be true 
along the Y-axis too, but the actual gridding Intervals along the Y-axis 
close to the poles are smaller than those close to the Equator. It Is 
obvious that this data base structure does not have equal -area and evenly 
distributed elements. 

The simplicity of the coordinate system allows such a data base to 
be capable of quick and direct-retrieval of data with minimal computa- 
tions. All data In the data base may be Indexed by directly calculated 
addresses. The detenril nation of data base addresses Is straightforward. 
Once the latitude-longitude coordinates are available, the transformation 
to the rectangular storage coordinates Is Just a matter of relabelling 
the coordinates. The binary serial address can then be obtained by 
simple table look-up procedures. 

2.1 Telescoping Addressing Scheme 

The specific addressing scheme discussed In this section Is based on 
the binary division scheme. The Earth's surface Is divided Into a number 
of square regions. Each region Is further divided, to the requisite 
resolution level, by a two-dimensional binary grid, as shown In Figure 
2-2. On each level of division, the areas are divided Into quadrants, 
which are labeled by a 2-b1t binary number. Each level of division Is 
Indicated by the addition of two binary bits to the least-significant end 
of the address. This binary address Is the serial location of a point In 
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the square region. An address to the third-level of division Is pre- 
sented In Figure 2-3 as an example. The expression of an address as a 
single bit string allows storage of addresses as single machine words and 
the expandiblllty and generality of the serial string permits the use at 
any resolution level by "telescoping". 

The surface flags are organized In area blocks within area blocks In 
a telescopic fashion. The main purpose Is to reduce the total data vol- 
ume required. This Is done by examining the homogeneity of large area 
blocks. For example, a large area block over land only requires the 
storage of a single surface flag. In one experiment, the ntxnber of homo- 
geneous area blocks over land was determined for different block sizes. 

A 25 % reduction In data volume was shown to be possible. Figure 2-4 
shows the plot of the heuristic data. 

Telescoping procedure can be thought of as a tree-like structure. 

The procedure Involves a sequence of partitioning area blocks Into finer 
resolution. And. It Is characterl zed by the fact that each partitioning 
step Is subjected to examinations and decisions before further partition- 
ing continues. When the structure Is diagrammed to show the h1e*!*archy 
among the procedure. It exhibits a characteristic tree-like aspect. 

Figure 2-5 Is an example to show the conceptual telescoping pro- 
cess. In this example, the address of a resolution element Is broken 
down Into four parts. The first part of the address. Ai . points to a 
large area block labelled AREAl. The bit string ^2 addresses a smaller 
area block AREA2 In which AREA2 Is contained In AREAl. Similarly. As 
and Ai» address AREAS and CELL, respectively. In actual storage. AREAS 
might be structured as a logical record containing resolution cells. At 
each of the telescoping steps, a decision Is made based on the homogen- 
eity of the area examined. The process continues If the area block Is 
heterogeneous. This corresponds to the continuous processing of the tree 
at a lower level. The addressing process ends If the area block being 
examined Is homogeneous or the tree level with the desired resolution Is 
reached. 
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Figure 2-5 

The Conceptual Data Retrieval Procedure by Telescoping 

Four levels are used in this example to address a 
resolution cell of desired si 2 e. The tree-like 
structure representing the process is also shown. 






Bi'sinkss and Tecmnouhucal Systems, isc 


2.2 Cost Analysis of Telescoping 

There are two parameters to be determined for the telescoping 
retrieval procedure. The first parameter Is the number of telescoping 
steps required to achieve 7 km resolution and the second Is the area 
block size at each step. An optimization procedure Is formulated to 
determine these parameters. It Involves the formulation of cost functions 
at each and every telescoping level and the analysis of some empirical 
data. 


Telescoping Is based on the fact that there exists large homogeneous 
area blocks, that Is, blocks consist of a single type of surface. It can 
be easily seen that the percentage of heterogeneous area blocks on the 
Earth's surface decreases with the area block size. Heterogeneous area 
blocks are area blocks that contain different surface types (for example, 
blocks contain land-ocean boundary). The relationship between the amount 
of heterogeneous area blocks on the Earth's surface and different block 
sizes is empirically determined. Figure 2-6 plots the data. 

The x-axIs, , represents the number of partitions, 7 km apart, on 
a side of the area block and the y-axis, f(M) , represents the percen- 
tage of heterogeneous area blocks of the total Earth's surface area. The 
plot reveals a linear relationship between M and f(M) when M<400 . 

As M approaches zero, the curve Intersects the y-axis at f(M) » 0.1 . 
In other words, there are about 10% heterogeneous cells (I.e., land-ocean 
boundaries) on the Earth's surface when the resolution approaches 7 km. 
This linear relationship Is Incorporated Into the optimization anclysis 
described as follows. 

Let n. be the number of partitions on a side of a square area 
*3 ^ u 

block at the j telescoping level, and N be the total number of 7 km 
partitions along the latitutde, or half ^he perimeter of the Earth. N 
can be written as: 
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Figure 2-6 M vs. f(M) 
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1 

N » n, n, ••• n. ■ n n. 

^ j-1 

where 1 Is the required number of telascoping levels to achieve 7 km 
resolution. N Is a constant expressed as: 

N - (y) 2ffR/7km 

• 28b9.75 


where R » 6371 km, the Earth's radius. 

Further define a function f describing the plot In Figure 2-6 by 

f(x) ■ a + b X 

where a Is the Intercept and b Is the slope. 

The fractional amount of heterogeneous area blocks encountered at 

f’h 

the k telescoping levels can be expressed as: 


f(“) = a + 

m. 


b(^) 


k 

where m. = n n. . 

^ j-1 ^ 

The cost to telescope 1 levels Is then expressed as: 

C = cost of level 1 + cost of level 2 ••• + cost of level 1 
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c • "1 * "1 "2 "2 - "? 


1 ^ O 1 M 

C - i [ \ n/] [ n f(^)] . 


, I '■ Mi J L " 

-1 j-1 J k-I 


k 


To minimize the above cost function, the extrema have to be 
located. This involves the solving of the following set of equations for 

"j • 

5r 

■ 0 where j*l,***i 


To illustrate the process, a few simple cases of evaluating the 
nj's are detailed below- 

Case 1 ; i*l , implies no telescoping is applied. 


2 

N ■ n^ , and C » n^ 

nj ■ 2860 

C - 8.178 X 10® 


Case 2 : 1*2. Two levels of addressing 


N = nj n^ I and 


.2 . „2 2 -/N . 

"i * "i "i 


„2 . „2 „2 

”l ”l ”l 
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To solve: - o • 2n. - b5!- 

1 fij 


'1 


The results; Oj - N(b/2)^/^ - 2C9 


nj - ( 2 /b )^/3 - 10.63 


C - 1.103 X 10 ® 


1*3. Three levels of addressing 
N ■ nj 02 03 , and 


. nj . ,2 ^ ^2 „2 „2 

■ "1 * * "1 n| n| (a*b^) (a* 


To solve: 4^ » 0 


«C 


= 0 


<?c 


The results: - [5^' . 3 bV j'/’ 3^^ 

^ (b/ 2 )^'^ 


11 


”2 “ nj (b/ 2 )^/^ 


1.99 


2-n 
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1 


"3 




in.63 


C 


- 0.17x10® 


Case 4 ; 1 ■ 4 . Four levels of addressing 

N “ P 2 n^ , and 


■ "l * "l "2 * "1 "2 "3 * 


2 .2 2 ,/N 


nf n\ nl f(-^) f(^) f(^) 
1 2 3 4 ' nij ' ' m2 ' m^ ' 


The results: n^ ■ 72.25 

n2 ■ 1»87 
nj - 1.99 
n^ ■ 10.63 

C - 0.034x10® 

To summarize, the following table lists results of the first four 
cases. Although the n^'s are optimally computed, they are re-assigned 
to Integers which are powers of two. This enables the binary serial 
address to be used most efficiently. In all cases, the number of parti- 
tions at the last telescoping level Is 12. The area block Involved In 
the last level of partition can be considered as a logical record. In 
this case, the record size Is 12 x 12 » 144. 
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i 

Computed n^ 

Reassigned n^ 

Cost 

1 

2857.7 

- 

8.178 X 10® 

2 

269.0 10.6 

256 12 

1.185 X 10® 

3 

135.1 1.9 10.6 

128 2 12 

0.179 X 10® 

4 

72.3 1.9 2.0 10.6 

64 2 2 12 

0.035 X 10® 


The cost of the five cases is Illustrated graphically in Figure ?-7. 
Case 4 is recommended. It consists of four telescoping steps. At the 
first level, the World Surface Map is partitioned into ? x 64 x 64 ■ 8192 
square areas employing 13 bits of address. The second step partitions 
each and every one of the 8192 areas into four squares, adding another 
two bits to the address and resulting in a total of 4 x 8192 ■ 32,768 
area blocks. The third telescoping level is the same as the second one. 
The last step addresses the 144 elements within a record and requires an 
address of eight bits. Figure 2-8 summarizes this recommended telescop- 
ing procedure. 


T 

</) 

O 

o 


1 



Figure 2-7 The Cost of Telescoping 
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3.0 EXAMPLE OF A DATA BASE STRUCTURE BASED ON THE SPHERICAL ICOSAHEDRON 


3.1 Introduction 


The spherical Icosahedron data base provides some useful properties 
which are valuable In many applications. In this model, the sphere Is 
visualized as a spherical Icosahedron as Illustrated In Figure 3-1. This 
spherical Icosahedron Is obtained by radially projecting the edges of an 
Inscribed Icosahedron. 



Figure 3-1 Spherical Icosahedron 

3.2 Geometric Properties of a Spherical Icosahedron 

A regular n-polyhedron Is a geometric solid whose faces are regular 
n-sided polygon. It Is well known that there can be only five regular 
polyhedrons as shown In Figure 3-2. 


3-1 




BusiyEss 4M) TECwmuMiHUL Systems, Inc. 



Regular Polyhedrons 

Number and Type of Faces 

Tetrahedron 

4 equilateral triangles 

Octahedron 

8 equilateral triangles 

Icosahedron 

20 equilateral triangles 

Hexahedron (cube) 

6 squares 

Dodecahedron 

12 regular pentagons 


Figure 3-2 The Five Regular Polyhedrons 


Each regular polyhedron can be Inscribed In a sphere, and the edges can 
be radially projected onto the sphere forming "spherical polyhedrons". 

The spherical icosahedron formed out of 20 spherical triangles Is the 
division of the sphere having the greatest number of regular pieces. For 
the triangles of the spherical Icosahedron, Illustrated In Figure 3-3, 
the vertex angles have the value a ■ 2^/H . 



Figure 3-3 Face of Spherical Icosahedron 







Biyii\K\s 4.Mi TKiMyouHin.41. SrsTKMS, l/vr: 


It can be shown that the sid'^’ Lp of ^ plane equilateral triangle 
of the regular Icosahedron Is given by 

■ 2R sin 4 
P < 

where R Is the radius of the sphere and “a" Is the side of the spher- 
ical triangle given by: 


a 

a - 2 cos- or 


Lj ■ aR . 


The area of the spherical triangle, A , and of the plane equilateral 
triangle, Ap , can be canputed by: 

■ (^) 4«R^ , and 


A„ ■ ~T • 
P 4 p 


Using 6371 kilometers as the radius of the Earth, the following list 
gives useful values Involving the above equations. 


R 

a 


6371 km 

1.2566 

1.1071 

6698.8659 km 
7053.6527 km 
19431360.3 km^ 
25503223.6 km^ 
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3.3 Spherical Icosahedron Data Base 

In constructing a data base structured about the Icosahedron solid, 
the process Involves the division of the 20 spherical triangles Into 
smaller triangles of nearly equal area. Consider the partitioning of the 
plane triangle of the regular Icosahedron Into equal-area equilateral 
triangles by successive processes of bisection of the sides. After four 
such bisection processes, the original triangles are partitioned Into 64 
smaller triangles as shown In Figure 3-4. Fach additional bisection 
quadruples the number of partitions. The data base element can be 
defined In form of a hexagon and each hexagonal resolution element has 
twice the area of a small triangle. 



Figure 3-4 An Example of Partitioning an Icosahedron Triangle 
(Resolution elements are defined as hexagons.) 
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Using the -bove example shown 1n Figure 3-4, a basic triangle of 
7053.65 kilometers on all sides consists of 32 resolution elements at the 
resolution of about 881.7 kilometers In width and 7C5975.7 square kilo- 
meters In area. To approach the resolution of 7 kilometers, the funda- 
mental length of the basic triangle has to be partitioned Into 1024 
parts, bach hexagonal resolution element corresponds to about 48.6 
square kilometers. Figure 3-5 shows other alternatives possible for 
resolution element sizes. 

From Figure 3-4, It Is obvious that equal-area elements on the plane 
triangle do not radially project as equal-area elements on the spherical 
triangle. For example, those elements near the center of the plane tri- 
angle have larger projections than those elements near the edges and ver- 
tices of the plane triangle. Hence, If a triangular grid of equal -area 
elements Is first constructed on the plane triangle. It Is then necessary 
to distort this grid Into a curvilinear network so that the elements near 
the center are smaller than those near the edges. The distortion Is such 
that when the curvilinear elements are projected radially, equal-area 
elements are again obtained on the spherical surface. The desired 
sequence of transformation Is Illustrated In Figure 3-6. The mathemat- 
ical details of deriving these transformations are discussed In the fol- 
lowing section. 

3.4 Mathematical Formulation of the Mapping Function 

First, consider a plane surface subtended by a spherical surface 
with radius R . Let be the vector from the center of the sphere to 

the given plane. As shown In Figure 3-7, let dAp be an area element on 
this plane, and let 7 be the vector from the center of the sphere to 
the area element dAp . Let dA^ be the spherical area element obtained 
by projected dAp radially onto the sphere. Then, It is not too diffi- 
cult to show that the area element dA^ on the sphere Is related to the 
area element dAp of the plane triangle via the relation; 
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SIDE 

PARTITIONS 

NUMBER OF 
HEXAGONS 
(million) 

WIDTH OF 
HEXAGONS 

AREA OF THE 
HEXAGONAL 
ELEMENT 

25o 

0.65 

27.6 km 

778.3 km^ 

512 

2.62 

13.8 km 

194.6 

1024 

10.48 

6.9 km 

48.6 

2048 

41.94 

3.4 km 

12.2 

4096 

167.77 

1.7 km 

3.0 


Notes: 

• Side partition is t’he number of partitions of a side of a basic 
icosahedron triangle. 

• Number of Hexagons is the total number of resolution cells on the 
Earth-oriented spherical surface and is computed by 1/2 (side 
partitions)^ x 20. 

• The width of a hexagon is measured as the shortest distance 
between two parallel sides. 

• The area of a hexagonal element is computed by (Total Earth's 
Surface Area) (Number of Hexagons) 


Figure 3-5 Various Resolution Sizes of the Icosahedron Data Structure 
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Plain Triangles with Equal Areas by Partitioning 



Plain Triangles with Distorted Grids by Mapping 



Spherical Triangles with Equal Areas by Projecting 

Figure 3-6 Sequence of Transformation 
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Figure 3-7 Relation Between Plane and Spherical Area Elements 
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dA 


s 


cos^ f”.F ) 



where f^»rQ) denotes the angle between r and . 


This relation Is used to derive the mapping function In transforming 
to equal-area elements on the sphere. The mapping function Is expressed 
In a new coordinate system (u,v,w) defined as lines of constant Inter- 
val forming the partitioning of the larger triangle (Figure 3-8). In 
particular, u-0 , v«0 , w»0 are lines forming the boundary of the 
triangle; u»l , v»l , w*l correspond to the three vertices. Let 7 be 
a simple vector from the center of the Icosahedron to any point on the 
triangle and It can be expressed as: 

7 ■ uX + vF + wT 

where X , ¥ , and ¥ are vectors from the center to the three vertices 
of the triangle containing 7 . 7 can also be expressed In terms of the 
mapping function g(u,p) 


r 



9 


llw4lC 

W.P) 


where 

g(u,p) * I + Y(u-^) -3(1 -y) (u-y)^ 
- “I (1 -t)p + ^ u(u-y)^ 

p . ( u ^)2 + ( v -^)2 + ( w ^)2 


and 


Y^ “ 0.8290 
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Figure 3-8 Coordinate System in a Plane Triangle 
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Using this transformation, a data base can be constructed by pro- 
jecting, onto the sphere, the vertices of triangles obtained by succes- 
sive bisection of the faces of the icosahedron. These points are centers 
of nearly equal-area hexagon-like resolution elements. The I? vertices 
of the icosahedron are centers of pentagons. 

3.5 Organization of Data Base and Addressing 

The data base system adopted from the near equal -area division of 
the spherical icosahedron is most likely to be organized into a binary 
division scheme. The major objective is to provide "telescoping", that 
is, each finer level of resolution is contained in the previous resolu- 
tion. The addressing scheme being binarily indexable is capable of 
accessing each finer level of resolution. 

The mapping of the data base onto storage array may be accomplished 
in a number of ways. The classic method is mapping by columns and rows, 
as is done in the usual mathematical matrix notation. This column-row 
mapping may cause problems in the spherical icosahedron data base because 
of the nature of the data base structure - triangular structure. Other 
mapping techniques include scatter storage with one-to-one mapping and 
various consistent logical stringing techniques. 

■’’ht' specific scheme studied in this report is based on binary divi- 
sioii i leme. This mapping starts at the level of the faces of the spher- 
ical icosahedron, numbering these faces 1 through 20 as in Figure 3-9. 
Each face is divided, to finer resolution levels, by a triangular grid, 
as shown in Figure 3-10. On each level of division, the areas are 
divided into quadrants, which are labeled by a 2-bit binary number. Each 
level of division, N , is indicated by the addition of two binary bits 

to the least-significant end of a 5+N bit binary number. This binary 

2N 

address is the serial location of a triangular area in the 20 x 2 
elements location array. An address to the third level of division would 
appear in Figure 3-11. Figure 3-12 provides the possible levels of divi- 
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Figure 3-9 Icosahedron Face Numbering Scheme 
(A 5-bit string is required, from 00000 to 10011) 
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Figure 3-10 Binary Division Scheme 
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Figure 3-11 Example of Addressing by Binary Bits 
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LEVEL OF 
DIVISIONS 

SIDE 

PARTITIONS 

LENGTH OF 
THE SMALLEST 
TRIANGLE 

NUMBER OF 
ADDRESSING 
BITS REQUIRED 

1 

2 

3527.0 

5 + 2 

2 

4 

1763.5 

5 + 4 

3 

8 

881.8 

5 + 6 

4 

16 

440.9 

5 + 8 

5 

32 

220.4 

5 + 10 

6 

64 

110.2 

5+12 

7 

128 

55.1 

5 + 14 

8 

256 

27.6 

5+16 

9 

512 

13.8 

5+18 

10 

1024 

6.9 

5 + 20 

11 

2048 

3.4 

5 + 22 

12 

4096 

1.7 

5 + 24 


Note: The five bits are for addressing the 20 icosahedron triangles, 

from 00000 to 10011. The rest of the address (N bits) is for 
2N 

addressing the 2 triangular elements in the basic icosahedron 
triangle. 


Figure 3-12 Level of Division vs. Addressing Requirement 
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slons, their corresponding sizes (lengths of the triangles) and their 
corresponding address requirements. 

In developing the method of ordering data In the data base, there 
are obvious advantages of using a serial addressing scheme rather than a 
normal two-dimensional addressing scheme. The expression of addresses as 
a single bit string allows storage of addresses as single machine words, 
whereas a two-dimensional addressing scheme would require two or three 
words, Including one for the face number. Also, the expandibll Ity and 
generality of the serial string permits the use at any resolution level 
without the recalculation of the serial Index. 

The serial addressing string can be constructed simply by consider- 
ing the high order m olts as the record number, and the lower (N-M) 
bits as the address within record. The record number points to the tri- 
angular area containing the desired resolution element and address of the 
hexagonal cell Is represented by the lower bit string. A specific exam- 
ple of addressing a 7 km resolution cell Is provided below. Figure 3-13 
Illustrates the addressing procedure of a hexagonal resolution element 
and Its corresponding address. The first five bits of the serial string 
Is the face number indexing one of the 20 Icosahedron triangles. The 
next 12 bits represent one of the 4096 smaller triangles after six levels 
of division. The small triangle has a side of 110.2 km and Is considered 
as a record. The last seven bits Is the address of one of the 128 reso- 
lution cells within the record In a predefined serial fashion. One way 
to number the resolution cells Is Illustrated In Figure 3-13. The hexa- 
gonal resolution cell has a width of approximately 7 km. Two of the 
corner resolution cells and two of the elements are stored In an adjacent 
record. As shown In Figure 3-13, only half of the remaining side resolu- 
tion elements are associated with and stored In the logical record. 

Thus, with this scheme, there Is no duplication of storage. Other possi- 
bilities exist for the Indexing of hexagonal resolution elements forming 
logical records. There are 2(4'^-l)poss1ble alternatives are where n Is 
the number of elements. 
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Logical records consisting of 8, 3?, 128, 512, 2048, and 81<J2 are 
possible. Figure 3-14 shows relationships between various record sizes 
and their address structures at 7 km resolution. The basic concept of 
this bit string sequential structure Is based upon blocked data contain- 
ing resolution elements of a single record (triangular area). The effi- 
ciency of data retrieval Is depended on the chosen size of the block 
data. The size should be determined such that ntmber of storage accesses 
Is minimized. In one of the cases as listed In Figure 3-14, a triangular 
area (record) of 110.2 km on each side consists of 8192 resolution ele- 
ments at the resolution of about 7 km. If a sector of a disk track con- 
sists of 4096 data elements, only two accesses are required to retrieve 
data for the entire area (record). Storage based on the next larger 
record size, 32768 elements, would require at least eight access opera- 
tions. Hence, there Is great Interest In studying trade-offs between the 
record size and the retrieval efficiency. 

3.6 Summary and Future Research 

The properties of the spherical Icosahedron data base are summarized 
as follows: 

• Near-equal area cells 

• Even distribution on the Earth's surface 

• Variable resolution level by telescoping 

• Flexible addressing scheme 

• Minimal distortion 

• Flexible orientation of data base relative to Earth 

• Relatively simple filling and retrieving data. 

It Is clear that the spherical Icosahedron data base has some prom- 
ising and unique properties which ae valuable In meeting requirements of 
the Earth-referenced data base. An obvious next step Is to study method- 
ologies of rapid retrieval of data. This Is particularly Important for 
real time processing requirements. It Is also necessary to provide a 
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Figure 3-14 Relationships Between Record Sizes and 
' Address Structures at 7 Kilometers Resolution 
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means of rapid transformation between the geographic coordinates and data 
base structural coordinaes for addressing purposes. Also, demonstrations 
win be In order to validate all the proposed procedures Including the 
retrieval process, the filling method, and the serially transferring pro- 
cess of data to a "rectangular" configurated storage memory. 
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